Confined states and direction -dependent transmission in graphene quantum wells 
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We report the existence of confined massless fermion states in a graphene quantum well (QW) 
by means of analytical and numerical calculations. These states show an unusual quasi-linear de- 
pendence on the momentum parallel to the QW: their number depends on the wavevector and 
is constrained by electron-hole conversion in the barrier regions. An essential difference with non- 
relativistic electron states is a mixing between free and confined states at the edges of the free-particle 
continua, demonstrated by the direction-dependent resonant transmission across a potential well. 

PACS numbers: 71.10.Pm, 73.21.-b, 81.05.Uw 
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Recent studies have demonstrated the production of 
stable, ultrapure, two-dimensional (2D) carbon crystals, 
also known as graphene 0, 0, . These 2D crystals pos- 
sess unusual properties, such as unconventional quantum 
Hall effect 4, 5, 6, 7J and a strong electric-field effect 8j. 
A large part of these new properties are a consequence of 
the linear (in wavector) energy spectrum near the Fermi 
energy and are expected to lead to a new class of carbon- 
or graphene-based nanoelectronic devices. Previous the- 
oretical studies of relativistic fermions interacting with 
strong fields have indicated that the quantum behavior 
of the particles may differ considerably from the non- 
relativistic case 9j. In this paper we investigate the na- 
ture of electron states in graphene QWs and their quan- 
tized spectrum. 

Graphene layers consist of a honeycomb lattice of 
covalent-bond carbon atoms, which can be treated as two 
interpenetrating triangular sublattices, labelled A and 
B, and are often discussed in terms of unrolled, single- 
wall carbon nanotubes. The low-energy band structure 
of graphene is gapless and the corresponding electronic 
states are found near two cones located at unequivalent 
corners of the Brillouin zone . The low-energy carrier 
dynamics is equivalent to that of a 2D gas of massless 
charged fermions. Their behavior is governed by the 2D 
Dirac Hamiltonian [ill. H^. 
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where the pseudospin matrix a has components given by 
Pauli's matrices, p = {px^Py) is the momentum operator, 
and vf the effective speed of light of the system, which in 
this case corresponds to the Fermi velocity vp ^ 1 x 10^ 
m/s. The Hamiltonian (1) acts on the states represented 
by the two-component spinors ^ = [i/ja , '^b]^, where ipA 
and iIjb represent the envelope functions associated with 
the probability amplitudes at the respective sublattice 
sites of the honeycomb graphene structure. The low- 
energy spectrum of free carriers is ^ = ±hvF{k'^-\-kyy^'^ ^ 
with kx and ky the wavevector along the x and y axes, 
in the vicinity of the cones at the Brillouin zone; the 



+ (— ) sign refers to electron (hole) bands. Equation (1) 
also implies that the carriers are chiral particles, with the 
pseudospin aligned parallel (antiparallel) to the direction 
of propagation of the electrons (holes). 

Representing the effect of an external electrostatic field 
by an external potential U and including a diagonal ef- 
fective mass-like term mvp leads to the Dirac equation 



[vF{a-p) 
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The term oc mvjp creates a gap in the dispersion and 
may arise from spin-orbit interaction or from the cou- 
pling between the graphene layer and the substrate [l^ . 
For a circularly symmetric potential with m = 0, the 
solutions inside the potential well match free-particle so- 
lutions outside, therefore ruling out bound states 
This is caused by the conservation of the chirality in the 
interaction with the potential and the absence of a gap in 
the spectrum and can be understood as a manifestation 
of a relativistic tunneling effect first discussed by Klein 
I15L for one-dimensional (ID) potentials, in which 
Dirac fermions can propagate to hole states across a steep 
potential barrier without any damping. For massless par- 
ticles this tunnelling is expected to occur for any value of 
Uq. However, as we show below, for a ID potential a finite 
value of the momentum parallel to the potential barrier 
can suppress this tunnelling and thus allow the confine- 
ment of electrons. Very recent studies have demonstrated 
the confinement of electrons in a graphene strip 17]. In 
this case, in order to obtain the confinement the authors 
assumed a position-dependent effective mass for the par- 
ticles. This assumption does not permit the observation 
of Klein tunnelling and of the momentum-dependent re- 
flection and transmission. Therefore, the confinement 
in this case is qualitatively different from ours specified 
below. In order to demonstrate the confinement in an 
electrostatic quantum well, we consider a zero or con- 
stant effective mass throughout the system and first a 
ID square-well potential U{x) = Uq 0{\x\ - L/2), Uq > 0, 
cf. Fig. 1, which allows an analytical solution for the 
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FIG. 1: A square quantum well on a graphene layer. 



eigenstates and sheds light on some general features of 
the problem. Later on, we consider a parabolic confine- 
ment. 

With momentum conservation in the y direction, we 
look for solutions in the form iljc{x,y) = (j)c{x)e'^^yy , 
C=A, B, and obtain 

d(t)B/d^^ p(t)B = i{e - u - A)(l)A, (3) 

d(t)A/d^- P(l)A = i(e-ii + A)0B, (4) 

where ^ = x/L, /3 = kyL^ e = E L/hvp^ u = U{x) L/hvp 
and A = mvFL/h (for graphene hvp = 0.539 eV nm). 
Decoupling Eqs. (3) and (4) gives for (pA the result 



d^^ (e- 



A) 
+ A) 



(e - uf + A^ = 0, (5) 



where u' is the derivative of the potential. For a square 
well, these derivatives are Dirac delta functions. 

The character of the solutions depends on the value 
of which determines the sign of the last term on the 
left side of Eq. (5). The solutions are of three types: (i) 
traveling waves, which describe free electrons, free holes, 
as well as mixed states that occur due to the Klein tun- 
neling of electrons to holes outside the potential well; (ii) 
standing waves, which for massless fermions arise only 
from finite values of 13 above an energy-dependent cut- 
off and decay exponentially in the barrier regions; (iii) 
tunneling waves, which are oscillatory outside the well 
whereas inside it they are combinations of exponentials 
with real exponents; these correspond to holes that un- 
dergo ordinary tunneling across the potential well. Type 
(ii) solutions occur in energy and wavevector ranges for 
which there are no hole states available at the barrier 
regions. This suppresses the Klein tunneling, since it de- 
pends on the electron- hole conversion at the interface. 

In this work we focus on type (ii) solutions which de- 
scribe electron states confined across the well and prop- 
agating along it. Their energies are in the region delim- 
ited by the curves E = [{hypkyY + rin?VpY^'^ + and 



E = [{hypky)'^ -\-m'^Vp]^^'^ . At smaller wavevectors, tun- 
neling across the barriers introduces a cut-off in the spec- 
trum for £^ < —[{hypky)'^ -^m'^Vp]^^'^ -^Uq. For confined- 
states, the spinor components decay exponentially in the 
region ^ < —1/2. Then the A component can be written 
as (/)a(0 = ^1 Substituting 0a (0 Eq. (4) we find 
0s (0 = ^/-^i e"^, with /_ = {p-a)/{e-uo^A) and 
the decay constant a given by a = — (6—1x0)^+ A^]-^/^, 
where uq = UoL/hvp- 

The solutions (j)A and 0b for |^| < 1/2 are of the type 

0a(O = C2 cos«) +1^2 sin«), (6) 
<I^b{0 = [V(e + A)] C2[/3 cos«) + sin«)] + 

D2[P sm{K^) - kcos{k^)]\, (7) 

with f^'^ = - - A^. For ^ > 1/2 the solutions 
are similar to those for ^ < —1/2 but with a nega- 
tive exponent: 0a (0 = ^3 e""^, 0b (0 = ^/+ ^3 e""^^, 
where /+ = {/3 -\- a) / {e — uq -\- A) . It should be stressed 
that, in contrast with the non-relativistic case, the spinor 
components are neither even nor odd functions, despite 
the symmetry of the potential. This symmetry, how- 
ever, is refiected in the probability density p = ^'''^ = 
0a(OVa(O + 0b(OVs(O .14] , which is an even func- 
tion. Moreover, for a step potential the derivatives of the 
spinor components are not continuous because u' in Eq. 
(5) is a delta function. This can be demonstrated by con- 
sidering the continuity of the y component of the proba- 
bility current, jy = 'U/^V^y^, across the potential inter- 
face: using Eqs. (3) and (4) we obtain {u-^ = uo/{e-\-A)) 

^AA'^m = (1 - M+)(/''^^(l/2) + m+/3(/.a(1/2), (8) 

where the arrows indicate the limiting values from the 
left and right of the interface. Notice in Eq. (8) that, 
even for large values of A, a continuous derivative of 0a 
may be assumed only for uo(3 << A. 

Requiring the continuity of 0a and 0b at ^ = —1/2 and 
^ = 1/2 we obtain the following transcendental equation 
for the energy eigenvalues 

S- (e, A +1) 5+(e, A +1) + S- (e, A -1) ^+(e, A -1) = 0, 

(9) 

where S±{e,(3,s) = - /±(e + A) - si^S^' and S = 
tan(A>:/2). The non-relativistic limit can be obtained us- 
ing e = ec + A, where Cc corresponds to the classical 
energy and considering the limit A >> ec, to give 



/±(e + A)^(/3±a)(l + r), 



(10) 



where a ^ [/^^ + 2A(iio - ec)]^/^ 1^ ~ (2Aec -/3^)^/^ and 
r = uq/2A. Equation (9) then becomes (R = hz/2) 

[a(l + r)-A>:tan^] [a(l + r) + /^ cot ^] -2(/3r)^ = 0. (11) 

For r << 1 and /3uo << A we recover the familiar tran- 
scendental equation for a non-relativistic QW. In this 



limit a non-zero value for p is equivalent to a simple shift 
of the energy scale ^ — /5^/2A and the spectrum 
of the confined states becomes a set of nested parabolas. 
On the other hand, Eq. (11) shows that, even for mas- 
sive particles, the QW spectrum does depend on the y 
component of the momentum, in contrast with the non- 
relativistic results. Thus, a significant modification of 
the parabolic spectrum occurs as P increases. 

Equation (9) was solved numerically. The results are 
shown in Fig. 2 for [/q = 50 meV, L = 200 nm, and 
A = 0. The dashed lines delimit the continuum region, 
which corresponds to free-electrons {E > hky -h Uq) with 
energies greater than the barrier height, and free- holes 
{E < —hvFky-\-Uo) that propagate in the system by 
means of the Klein tunneling mechanism. The cut-off 
at low wavevectors thus arises due to the conversion of 
confined electrons to free holes. For large values of ky 
the dispersion branches are given approximately by 



E = hvF[{n7r/Lf^kl]^/^ 



(12) 



where n is an integer. For any given ky^ the accuracy of 
this approximation improves as L increases. The lower 
inset shows (a) (pA (solid curve) and (dashed curve) 
for the confined state, with ky = 0.03 nm~^, shown by 
the solid triangle and (b) the corresponding probability 
density in arbitrary units. The plot clearly indicates a 
discontinuity in the derivative of the spinor component 
functions at the barrier interfaces. The vertical dotted 
lines indicate the walls of the well. The upper inset 
shows the effect of the mass, with mvp = 10 meV. The 
dashed lines again represent the limits of the free-particle 
continua. In this case, confined states are allowed, for 

= 0, in the range uq — A < e < uq -\- A. This energy 
range broadens as ky increases and remains constant for 
ky > (ix§/4 — A^)^/^. At lower energies, there is again a 
cut-off, due to the Klein tunnelling at the barriers, which 
disappears for 2 A > i^o- 

Next, we consider a QW with a parabolic potential pro- 
file U{x) = Uo{2x/L)^ for \x\ < L/2 and U{x) = Uo for 
\x\ > L/2. Figure 3 shows the spectrum of confined states 
obtained from a numerical solution of Eqs. (3) and (4) 
for Uq = 50 meV and L = 200 nm. The results are qual- 
itatively similar to those of the previous case, but now 
with the eigenvalues being approximately equally spaced 
for large wavevectors. 

An essential difference with non-relativistic electrons, 
evident in all cases, is the appearance of new confined 
states at the edges of the continua, where the quan- 
tized electron branches intercept the free-particle regions. 
Thus, by an adiabatic increase in ky one can transform 
a free-electron or a free-hole state into a bound electron 
state. This occurs because the presence of the barri- 
ers allows a mixing of electron and hole states with the 
same energy and y component of momentum. As a re- 
sult there is constructive interference between confined 
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FIG. 2: Spectrum of confined states in a graphene square QW 
vs ky for Uq = 50 meV, L = 200 nm, and mvp — 0. The lower 
inset shows (a) and for the state shown by the solid 
triangle and (b) the related probability density. The upper 
inset shc2^§ the (^ffegt of, a nc/n-z^ro Tjias^, for; in,yjr = 10 meV. 
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FIG. 3: As in Fig. 2 but for a parabolic QW. 



states and unbound electron or hole states that are reso- 
nantly transmitted across the QW. We demonstrate this 
by calculating the transmission coefficient of electrons 
incident on a square well. Consider the propagating so- 
lutions ipA{x,y) = (t)A{x)e'^^yy , with 



(t>A{x) 



e < -1/2, 



A^e''^^ -h B2e-''^^ -1/2 < ^ < 1/2, (13) 



Ase 



^ > 1/2, 



where a = [{e — u^Y — P'^ — A^]^/^; the solutions for 0^ 
are obtained as in the previous calculation. Then, the 
transmission coefficient is obtained as T = |^3p, where 

, ^ (g+-fl-)(/+-/-)e-^" 

' {g+ - U){g. - /_)e2- - (5+ - f_){g_ - /+) ' 

(14) 

g± = {f3± ia)/{€ + A) and f± = {f3± m)/(e -uo + A). 
A {ky,a/L) contour plot of the transmission T is shown 
in Fig. 4 for Uo = 50 meV and L = 200 nm. As seen, 
T depends on the direction of propagation and displays 
an oscillatory behavior. As a — > 0, T reaches a series of 
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FIG. 4: Contour plot of the transmission coefficient of elec- 
trons incident on a graphene square well, with energy E > Uo, 
as a function of a, for Uo = 50 meV, L = 200 nm and m = 0. 



maxima for values of /3 that coincide with the wavevec- 
tors for which mixing occurs, cf. Fig. 2. Notice that for 
a significant range of incident angles T is always equal 
to 1. This includes the case of nearly normal incidence, 
ky ^ 0^ and is in sharp contrast with the non-relativistic 
case in which T exhibits periodic maxima equal to 1 as 
a function of kx. A similar direction-dependent trans- 
mission through graphene barriers was reported recently 
18]. A direction-dependent transmission is also possible 
for non-relativistic electrons tunneling through magnetic 
barriers jl9j . 

The y components of the momentum for which mixing 
is allowed correspond to confined states for which the 
asymptotic limit a ^ applies. This yields the condition 
sin (/^) = or = nTT, where n is an integer. Using the 
definition of and a gives 



\ 2uo 2 ) 



1/2 



(15) 



Since /^^ > 0, the values of n can be obtained from the 
condition ±(n^7r^/2iio — '^o/2) > A, where the + (— ) sign 
is associated with the upper (lower) continuum edges. 
From this condition we find that for Uq < 2mvp there is 
no mixing at lower energies, although it persists at the 
upper continuum edge and the minimum value of /3 for 
the mixing increases with A. 

A complementary way to see the direction dependence 
of the transmission T is shown in Fig. 5(a), with T plot- 
ted versus the angle of incidence = arctan (/c^/a), for 
different electron energies as indicated. The QW param- 
eters are L = 200 and Uq = 50 meV. Notice that for 
^ ~ 0, we have T ^ 1 in agreement with the ky ^ part 
of Fig. 4. In Fig. 5(b) we plot T versus the energy E 
for = 7r/3. As seen, T oscillates with the energy due 
to the resonance effect caused by the confined states (as 
in the Ramsauer-Townsend effect). The energies for the 
maxima of the transmission can be obtained from Eq. 
(15) as e = {rniY /2uq + uq/2. 

In summary, we showed that it is possible to confine 
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FIG. 5: (a) Transmission T vs angle of incidence 9 for differ- 
ent energies as indicated, (b) T vs energy E for — 7r/3. The 
maxima are marked by the numbers of the confined states. 



massless charge carriers by means of electrostatic poten- 
tials, due to the wavevector-dependent suppression of the 
electron- hole conversion at the potential barriers. We 
thus obtained the quantized spectrum of confined elec- 
tron states in graphene quantum wells (QW) as a func- 
tion of the y component of the wavevector. The results 
show a remarkable dependence of the eigenvalues on the 
momentum with a cut-off at low wave vectors. The rel- 
ativistic correction to the classical QW spectrum leads 
to a wavevector dependence of the number of confined 
states due to the electron-hole conversion at the contin- 
uum edges. Accordingly, such QWs must be treated as 
inherently 2D systems. This is further demonstrated by 
the directional dependence of the transmission shown in 
Figs. 4 and 5. Studying the resonance transmission of 
electrons across a QW with energies above the height of 
the confining walls, E > Uq, can probe the discrete levels 
which can be populated by tuning the Fermi energy of 
the system with the electric-field effect [Ij . 
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